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We obtain the equations of motion for a ferromagnetic Bose condensate of arbitrary spin in the 
long wavelength limit. We find that the magnetization of the condensate is described by a non- 
trivial modification of the Landau-Lifshitz equation, in which the magnetization is advected by the 
superfluid velocity. This hydrodynamic description, valid when the condensate wavefunction varies 
on scales much longer than either the density or spin healing lengths, is physically more transparent 
than the corresponding time-dependent Gross-Pitaevskii equation. We discuss the conservation 
laws of the theory and its application to the analysis of the stability of magnetic helices and Larmor 
precession. Precessional instabilities in particular provide a novel physical signature of dipolar 
forces. Finally, we discuss the anisotropic spin wave instability observed in the recent experiment 
of Vengalattore et. al. (Phys. Rev. Lett. 100, 170403, (2008)). 
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I. INTRODUCTION 

Perhaps the most dramatic way in which Bose con- 
densation in the alkah gases differs from its counterpart 
in '*He (or, for that matter, in conventional s-wave su- 
perconductors) is that the condensed particles have non- 
zero spin [T] . Magnetic trapping results in a gas in which 
the spin state of the atoms can be described using the 
adiabatic approximation, with the atoms remaining in a 
particular hyperfine level relative to the local magnetic 
field as they move around the trap. Although there are 
several interesting consequences of the non- uniformity of 
the field jSj, far richer behavior results when atoms are 
optically trapped, allowing the full consequences of rota- 
tional invariance to be realized. The experimental prepa- 
ration of a long-lived gas in a particular hyperfine mul- 
tiplet may be more difficult for some atoms than others, 
but the theorist is nevertheless called upon to answer the 
question: what are the properties of the higher spin Bose 
condensates, and how will they manifest themselves in 
the ultracold laboratory? 

Starting with Refs. |3j IH [5l [6], a number of investi- 
gations have explored the possible magnetic phases of 
these condensates. The dynamics of such an ordered 
phase is generally described by a non-linear equation for 
the motion of the order parameter on some manifold of 
symmetry- broken states, but this line of thought has not 
been much pursued. This is partly due to the existence of 
a dynamical description of Bose condensates valid in the 
dilute limit, namely the time-dependent Gross-Pitaevskii 
(GP) equation, which may be straightforwardly extended 
to the multicomponent case. A hydrodynamical descrip- 
tion of the 'slow' degrees of freedom - generally the order 
parameter and any conserved quantities - is nevertheless 
desirable both for physical transparency and simplicity, 
at the modest cost of eliminating certain uninteresting 
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high frequency motions. The goal of this work is to de- 
velop and apply such a description for the case of a fer- 
romagnetic condensate. 

The structure of the remainder of the paper is as fol- 
lows. In the next section we derive the equations of mo- 
tion valid in the long wavelength limit starting from the 
Gross-Pitaevskii Lagrangian for a ferromagnetic conden- 
sate. After discussing the conservation laws of the theory 
the equations of motion are used to study the stability of 
magnetic helices (Section III I, Larmor precession in the 
presence of dipolar forces (Section IV ) and spin waves 
in the presence of Larmor averaged dipolar forces (Sec- 
tion |v]) . The relevance of the last calculation to the ex- 
periment of Ref. [71 is briefly discussed. 



II. EQUATIONS OF MOTION 

Our starting point is the Gross-Pitaevskii Lagrangian 
density {h = m — 1) 
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where <i> is a 2s -I- 1 component spinor and S the spin- 
s matrices. The quartic terms are appropriate to the 
description of interatomic interactions in a spin-1 con- 
densate such as ®'^Rb, where C2 < favors ferromag- 
netism [3111]. The description of higher spin condensates 
requires more parameters but the phase diagram always 
includes a ferromagnetic phase [5j [6] . 

We are going to work in the low energy limit where 
both interaction terms are fully satisfied. This is appro- 
priate to the limit where the condensate wavefunction 
varies on scales much longer than either the density or 
spin healing length, or equivalently, the superfluid ve- 
locity is small compared to the speed of propagation 
of sound or spin waves. This is analogous to the in- 
compressible limit used to describe normal fluids at low 
local Mach number. The density interactions demand 
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p(r) = $t(j) = const., which we set equal to unity from 
now on, while the spin interactions, assumed ferromag- 
netic, demand that the polarization be maximal at each 
point, so that the spinor <I> is a spin coherent state. As- 
suming that the interactions dominate all other terms in 
the Hamiltonian allows other terms describing e.g. Zee- 
man and dipole-dipole interactions to be included in a 
controlled way by simply evaluating them on the con- 
strained ferromagnetic manifold. A similar approach to 
the static case was introduced in [3 [9] . 

In this limit the equations of motion for the unit vector 
n(r, t) describing the local magnetization sn — <i>^^S<i> will 
be shown to be 



Dn 1 . 

n X V^n = 

Dt 2 



(2a) 



V-v = 0, 



V X V = -ea/37f^aVr^^ x Vn-^, (2b) 



where D/ Dt denotes the usual Eulerian derivative dt+'v- 
V. The first equation is a modified Landau-Lifshitz equa- 
tion (LLE), that accounts for the advection of the mag- 
netization by the supcrfluid velocity v. The other two 
equations determine this flow from the condition of van- 
ishing divergence (incompressibility) and the Mermin-Ho 
relation [TD] that fixes the vorticity lj = V x v. The 
lines of vorticity coincide with the lines of constant n. 
Eqs. (2b I fix V up to a some potential contribution 
and V'^ip = 0- This should be chosen so that the nor- 
mal component of the velocity vanishes at the boundary 
of the flow. Note that the dynamics of the system de- 
pends crucially on the spin s, with the usual LLE being 
recovered in the s —^ limit. 

Variational principle The dynamical equations may be 
found by a variational principle through an appropriate 
parametrization of the constraint manifold. We use the 
parametrization 

where is a normalized eigenstate of n • S with eigen- 
value s. As expected there is some gauge freedom in how 
the overall phase of the state is apportioned between 
and the phase factor e'^. For instance, the velocity is 
given by V = Wd — a, where a = i<i>J,V<i>„. This vector 
potential depends on the gauge choice, though its curl 
does not 



V X a = iV$^ X V$ 



which is just the Memin-Ho relation. Incompressibility 
translates to the constraint 



V • (V6I - a) = 0, 



(3) 



while the Hamiltonian and Lagrangian Eq. Q take the 
form 



is(Vn)V^(V0-a)^ 



H = 

c = (at-e)-n. 



We obtain the equations of motion by a variation of 
the associated action, bearing in mind Eq. ([3|. The 
non-trivial part is finding the variation of the terms in- 
volving Op (/i — t,r), which depends upon n, without 
introducing a specific parametrization. This is accom- 
plished by writing the field strength dfj_ai, — d^a^ — 
—sEap^nad^ripd^n^ including a fictitious extra coordi- 
nate u with n(t, r, u) =n(f, r) -t- u5n{t, r). Then 



sn X df^n 



Sn + «9^<i>Jj(5$j, 



The variational derivative of the action with respect to 
is just the current = (1, v), so we have 

J drdt 5a^j^ = j drdt [sj^n x S^n • 5n - if^'^^S'^^d^.j^] , 

where integration by parts has been used to obtain the 
second term. The boundary term is equal to zero as Sn 
vanishes there. Since the current is conserved by the 
incompressibility condition, the second term drops out 
and we have our variation. A final subtlety is that when 
we vary Sn we also have to take into account the change 
S9 implied by Eq. ^ 

V^Se = V • (5a 

This variation is handled in precisely the same way, but 
now the boundary term is zero due to the vanishing of 
the normal component of the velocity at the boundary. 

Together with the variation of the first term of the 
Hamiltonian, the above variation readily gives the equa- 
tion of motion Eq. ( 2a I . 



Conservation laws Let us start by noting that Eq. (2a I 



can be understood as the conservation equation for the 
magnetization dtn + diJi ~ with the spin current 

Jj = nvi — -n X 9^11. 

Another conservation law follows by first considering the 
equation of motion for the vorticity a; = V x v [TT] 



Duji s 
~W " ~2 



Sijkdkdioij. 



Oij = -Sijdkn ■ dkn - d^n ■ djn. 

which can be used to check that the hydrodynamic im- 
pulse defined by 

I = - f r X uj, 

is a constant of the motion, as it is for both the LLE 
(where a; is defined from the Mermin-Ho relation) and 
the Euler equation for incompressible flow. If external 
forces act on the fluid, they give the rate of change of 
the impulse. For the case of rigid walls that we have 
considered the total momentum / v of the fluid naturally 
remains zero, with the reaction force of the walls being 



3 



transmitted instantaneously to the body of the fluid due 
to incompressibility |12j . 

Finally we have the topological helicity invariant 



n 



1 

167r2 



1 
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which again is known in both hydrodynamics and mag- 
netism [131 [13] . For a condensate it is equal to the in- 
teger Brouwer degree of the map g : M'' ^ SU (2) given 
by <&(r) = (7(r)<i>Q for some fiducial state ^q. Topological 
defects have been sought in the spin-1/2 case, but so far 
without success [131 [H] . 



III. INSTABILITY OF A MAGNETIC HELIX 

We now apply these equations to the analysis of the 
stability of helical configurations of the magnetization 

no(^) = cos 6 + sin 6 {e^ cos £, + Gy sin ^) , 

with = qz — uJot. The instabilities of these configura- 
tions were studied in the recent experiment Ref . [^ • Since 
there is only a single wavevector present a; = v = 0. Sub- 
stitution into Eq. (2a I gives ujq — ^q'^cosO. To analyze 



the stability of this configuration we introduce the or- 
thonormal Frenet-Serret frame {no, 61,62} consisting of 
axes parallel to no, d^no, and the vector perpendicular 
to both. These satisfy 











2] 














with curvature k = sin 6 and torsion r = cos 9. Small 
deviations from the helix are then written as n = no + 
77161 -|-7]262. It is instructive to first discuss the prediction 
of the usual LLE without the advective term, which gives 
the coupled equations of motion 

1 q2 
m = -2 ^^^2 - qcos 98^7]! - y sin^ 9r]2 

m = ^V^?7i - gcos6'92?72, 



with dispersion fi(A:) ~ q cos 9kz ± 
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revealing an instability for < fc < gsin6'. Note that 
the growth rate of the unstable modes is isotropic. The 
non-linear evolution of the helix in the one-dimensional 
integrable case is discussed in Ref. [TTI, but let us now 
return to the full equation of motion. At the linear level 
the vorticity is 

a;(r) = qs sin 9 (d^m^y - dyV2ex) 
u;(k) iqs sin 9 (kxBy - fcye^)ry2(k). 

Solving Eq. ( 2b I for the velocity gives 



v(k) 



qs sin 9 
fc2 



kyk^Gy 



^2(k), 



so that the linearized advection term is 

„2 „ ,,•2 



fc2 

(v • V) n = -g'^ssin"' 9-^ri2()i)ei, 
leading to the dispersion relation 



Q{k) — qcos 9kz ± 



fc2 



+ g2 sin^ ( 



kl 1 



Note that the growth rate of the unstable modes is now 
anisotropic, with the transverse modes always stable. 
This result has been checked in both the spin-1 and spin- 
1/2 cases by calculating the corresponding Bogoliubov 
modes of Eq. ([T]) before taking the incompressible limit 
cojC2 — > 00. An extensive analysis of the Bogoliubov 
modes of the helix configurations, not restricted to the 
incompressible limit, was recently given in Ref. |18j . 



IV. DIPOLAR INTERACTIONS AND THE 
LARMOR INSTABILITY 

We now turn to the question of the stability of Lar- 
mor precession. Absent dipolar interactions, rotational 
symmetry in the spin space means that the only effect 
of a magnetic field is precession at the Larmor frequency 
lul = gfisH. We will now show that accounting for dipo- 
lar forces renders Larmor precession unstable in general. 
The approximation of averaging the dipolar interactions 
over the Larmor trajectories, as was done in Ref. [T^ . 
for example, is therefore guaranteed to break down at 
sufficiently long times. 

The simplest geometry to consider is an infinite plane 
of thickness d. If the plane is perpendicular to the x-axis, 
the demagnetizing field for a uniform magnetization is M 
is — (M • Bx) Bx, which is on the order of 10^^ Gauss for 
a typical atomic gas. 

For a magnetic field in the z-direction, the magneto- 
static energy per particle is then scniag, with 

emag(n) = ^^^±nl + ujLriz, ujj_ = ^os{gfiBfp, 

where p is the density of the gas. The easy-plane 
anisotropy energy u}±_ due to the dipolar forces in on the 
order of /i x 10 Hz under typical experimental conditions. 
We may use the above energy as a local energy density 
as long as the magnetization varies on sufficiently large 
scales 3> d. As we will see in the next section, in 2d the 
long-ranged dipolar interactions have the non-analytic 
form kakf^/m, so that the leading deviation from 
this approximation is 0{kd) i.e. linear in the spin wave 
wavevector. Requiring that this part is small compared 
to the usual quadratic spin wave dispersion ujk = k'^/2 
then gives the two conditions uj±d ^ k ^ d^^. The 
equations of motion are then 



Dn 



- n X 



de 



mag 



On 



\v^n I = 



4 



Working in the canonical coordinates n^, (j), the preces- 
sion of a spatially constant magnetization no(i) obeys the 
equations 



ujl — Lo±_nz cos 



de 



magV 



1 - nl) sin 20 = - 



90 



(4) 



which may be solved exactly in terms of elliptic func- 
tions. For the case of small angle precession about the 
direction of the magnetic field, Kittel's classical result for 
the precession frequency is ^u>l(^ol + ppj . 

Introducing as before the Frenet-Serret frame we find 
the linearized equation of motion for the spin waves (since 
the zeroth order solution is constant in space the advec- 
tion term plays no role) |21j 

^ahdtrih = ^^^^70 + '^-L ('7a cos^ 0(t) - el{t)el{t)r]b) 

Where we have used the expression for the torsion T(t) = 
77.2 0, and the equation of motion Eq. Introducing 
the complex notation = (z^,Zfe), with Zk = -I- 
w/2fc)e'*®'='*'' we find 



[elel + 6^62] -t- w_L cos^ 



(5) 



Evidently there are resonances at ojl = ±Wfc . The size of 
the resonance region can be obtained with the trial solu- 
tion z = ae"^* -I- pe^C^'^k^'^^L-i^) ^ Solving the eigenvalue 
equation for lj for the lo^ ~ ojl resonance gives 
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Finite uj±^ opens regions of complex uj in the vicinity of 
the resonances where spin wave amplitudes grow expo- 
nentially. Generally we expect resonances at Wfc = ulol 
for integer n. The resonances are present for arbitrar- 
ily weak dipolar forces, and could provide a novel sig- 
nature of these interactions even in the limit uj^_ u>l ■ 
The long wavelength assumption used in the derivation 
is likely not necessary for the existence of an instabil- 
ity as both Eq. ^ for the uniform domain and the spin 
waves are common to the general case. The instability is 
a parametric resonance that results from their coupling. 
The comparatively large Zeeman fields used in current 
experiments do lead to the difficulty that the thermal 
spin wave amplitudes at wave vectors corresponding to 
energies ^ lol that provide the initial 'seed' for the insta- 
bility are very small. At present this is the main obstacle 
to the observation of the predicted effect. 



SPIN- WAVE INSTABILITIES WITH 
LARMOR AVERAGED DIPOLAR 
INTERACTIONS 



We illustrate the instability at lowest order in the dipolar 
interactions. At this order we evaluate the driving term 
in Eq. ([s]) fkit) on the unperturbed Larmor precession 



fk{t) 



[(1 + 772)6'"-* -(1-712)6 



If the parametric instabilities discussed in the previous 
section are unlikely to be seen in present experiments, it 
is natural to investigate the effect of the Larmor averaged 
dipole interaction. For a magnetic field in the 2;-direction 
this takes the form |19j 



5d 



^12 + 2/12 ^ 2z^2 



12 



n^(ri)77^(r2) - - (77-(ri)n-(r2) + n^(ri)77^(r2)) 



(6) 



where ri2 = ri — r2 and gd = fj-o (gsfis)^ ■ In this 
section we will show that the effect of this interaction 
is to generate a spin wave instability of a system with 
uniform transverse magnetization. The resulting spec- 
trum of unstable modes has a characteristic anisotropic 
structure in wavevector space (see Fig. [T|) similar to that 
observed in the recent experiment of Ref . ^ . This insta- 
bility has very recently been discussed from the Bogoli- 
ubov viewpoint in Ref. |22| . while we will focus on the 
long wavelength limit. This approach is not expected 
to be quantitatively valid everywhere in the unstable re- 



gion of wavevector space, but it captures the anisotropic 
structure near zero momentum exactly. 

Again we consider an infinite plane of thickness d in 
the y — z plane. To obtain the momentum space form of 
the dipolar interactions, it is convenient to make use of 
the formula 



'-12 + 2/12 2z^2 



47r 



' 12 



'ri2 



S^'\r,2), (7) 
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FIG. 1: Unstable mode spectrum Eq. (11 1 with uj_i_d^m/h'^ = 
0.25. The lower figure is a magnification of the region near 
the origin. 



with which we obtain 



H, 



dip 



Odd 

Stt J {2t:Y 



An 



(8) 



The second term in the square brackets arises from the S- 
function term of Eq. Q and accounts for the anisotropy 
energy introduced in the previous section: after Larmor 
averaging an easy-plane anisotropy an easy axis in the 
z-direction is obtained. The first term accounts for the 
finite wavevector correction to this picture and involves 
the function g(x) with asymptotes 



^ x> 1. 



(9) 



The low wavevector asymptote displays the non- 
analyticity characteristic of dipolar forces. 

From the point of view of energetics the easy axis 
term is the source of the instability that we will dis- 
cuss: a system that starts with transverse magnetization 
can lower its dipolar energy by ordering at long wave- 
lengths in the ±z directions. Incidentally, the asymptote 
Eq. ([9]) also shows that a tightly wound helical configu- 
ration with magnetization in the x ^ y plane and axis in 
the z-direction has the same dipolar energy as a uniform 
configuration aligned in the 2-direction, but of course the 
former has a lower kinetic energy. 

It is now straightforward to linearize the LLE about 
uniform transversely magnetized state (again we can ig- 
nore the advective term). Using the same notations as 
in Section |III| and considering only the case = of 
purely transverse magnetization, we obtain the coupled 
equations 



'?i(k) 



1 



klg{kd) 



^i(k) 



2 



UJ±(P 

47r 



klg{kd) 



^2(k), (10) 



which give for the unstable modes 




Ljj_\_d k1 



fc2 

y 



2 



LUj_d k1 
~^~k 



(11) 

where we have used the low wavevector asymptote of 
g{x). The imaginary part of the result Eq. ( [TT| is plotted 
in Fig.[l] The dimensionless dipolar energy uij_d'^m/ti^ — 
0.25 is appropriate to the experiment Ref. [7] with peak 
density po — 2.3 x lO-^'^cm"^ and d taken to be the 
Thomas- Fermi thickness 1.8 /im. 

The main qualitative feature of the spectrum of unsta- 
ble modes is the 'pinched' structure around zero wavevec- 
tor, with the boundary between stability and instability 



having the form k^ 



It seems plausible that the 



anisotropic power spectrum of the instability studied in 
Ref. [7j is a reflection of this structure. 

A remaining puzzle of that experiment concerns the 
observation that the anisotropic spin wave instability de- 
pends upon the creation of a helical twist of the trans- 
verse magnetization by a magnetic field gradient: no in- 
stability was observed in the case of a uniform system. 
On the other hand, the instability discussed above is 
hardly affected by the presence of the helix, not least 
because the characteristic scale of the the latter is much 
larger. It is then natural to suppose that the field gra- 
dients used to create the helix additionally generate spin 
wave fluctuations that provide a large initial 'seed' from 
which the instability grows. It would therefore be useful 
to study the same problem in a uniform system at higher 
temperatures, where such fluctuations may be present to 
begin with. 
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VI. CONCLUSION 

We have developed a long wavelength description of 
ferromagnetic condensates, which considerably simplifies 
the understanding of dynamical phenomena relative to 
the time-dependent GP theory. The utility of the ap- 



proach was illustrated with a number of examples that 
either have already been realized in the ultracold labora- 
tory, or may be in the near future. 

I would like to thank Mukund Vengalattore and Dan 
Stamper-Kurn for useful discussions of the experimental 
situation. 
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